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Abstract 

This paper describes topological searches for neutral scalar bosons produced in association 
with a Z'^ boson via the Bjorken process e+e~ S^Z^ at centre-of-mass energies of 91 GeV and 
183-209 GeV. These searches are based on studies of the recoil mass spectrum of Z" — > e"'"e~ 
and ^^/^~ events and on a search for S^Z*^ with — > and — > e"'"e~ or photons. They cover 
the decays of the into an arbitrary combination of hadrons, leptons, photons and invisible 
particles as well as the possibility that it might be stable. 

No indication for a signal is found in the data and upper limits on the cross section of the 
Bjorken process are calculated. Cross-section limits are given in terms of a scale factor k with 
respect to the Standard Model cross section for the Higgs-strahlung process e"'"e~ — > Hgj^Z'^. 

These results can be interpreted in general scenarios independently of the decay modes of 
the S". The examples considered here are the production of a single new scalar particle with a 
decay width smaller than the detector mass resolution, and for the first time, two scenarios with 
continuous mass distributions, due to a single very broad state or several states close in mass. 
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1 Introduction 



In this paper searches for new neutral scalar bosons S with the OPAL detector at LEP are 
described. The new bosons are assumed to be produced in association with a Z*^ boson via the 
Bjorken process e+e" S^Z^. Throughout this note, denotes, depending on the context, any 
new scalar neutral boson, the Standard Model Higgs boson Hgj^ or CP-even Higgs bosons hO m 
models that predict more than one Higgs boson. 

The analyses are topological searches and are based on studies of the recoil mass spectrum 
in — > e"'"e~ and events and on a search for S'^Z'^ events with — > e+e" or photons and 

Z° vD. They are sensitive to all decays of S*^ into an arbitrary combination of hadrons, leptons, 
photons and invisible particles, and to the case of a long-lived leaving the detector without 
interacting. The analyses are applied to LEP 1 Z^ on-peak data (115.4 pb~^ at ^/s = 91.2 GeV) 
and to 662.4 pb~^ of LEP 2 data collected at centre-of-mass energies in the range of 183 to 
209 GeV. In 1990 OPAL performed a decay-mode independent search for light Higgs bosons and 
new scalars using 6.8 pb~^ of data with centre-of-mass energies around the Z° pole Q. Assuming 
the Standard Model production cross section, a lower limit on the Higgs boson mass of 11.3 GeV 
was obtained. We have re-analysed the LEP 1 on-peak data in order to extend the sensitive 
region to signal masses up to 55 GeV. Including the data above the Z^ peak (LEP 2) enlarges 
the sensitivity up to mgo ~ 100 GeV. The S° mass range between 30 and 55 GeV is covered by 
both the LEP 1 and the LEP 2 analysis. 

The results are presented in terms of limits on the scaling factor /c, which relates the S'^Z'^ 
production cross section to the Standard Model (SM) cross section for the Higgs-strahlung pro- 
cess: 

o-gozo = A: • o-Ho^zo(mHO^ = mgo), (1) 

where we assume that k does not depend on the centre-of-mass energy for any given mass mgo. 
Since the analysis is insensitive to the decay mode of the S'^, these limits can be interpreted 
in any scenario beyond the Standard Model. Examples of such interpretations are listed in the 
following. 

• The most general case is to provide upper limits on the cross section or scaling factor k 
for a single new scalar boson independent of its couplings to other particles. We assume 
that the decay width is small compared to the detector mass resolution. In a more specific 
interpretation, assuming the S°Z° production cross section to be identical to the Standard 
Model Higgs boson one, the limit on k can be translated into a lower limit on the Higgs 
boson mas^. 

• For the first time we give limits not only for a single mass peak with small width, but 
also for a continuous distribution of the signal in a wide mass range. Such continua ap- 
pear in several recently proposed models, e. g. for a large number of unresolved Higgs 
bosons about equally spaced in mass ("Uniform Higgs scenario" ||3|]), or models with addi- 
tional SU(3)cxSU(2)lxU(1)y singlet fields which interact strongly with the Higgs boson 
("Stealthy Higgs scenario" 0]). These two models are described in more detail in the next 
section. 

2 Continuous Higgs scenarios 



^Dedicated searches for the Standard Model Higgs boson by the four LEP experiments, exploiting the predic- 
tion for its decay modes, have ruled out masses of up to 114.1 GeV fl. 
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2.1 The Uniform Higgs scenario 

This model, as described in Ref. ||3|, assumes a broad enhancement over the background expec- 
tation in the Mx mass distribution for the process e+e" Z^X. This enhancement is due to 
numerous additional neutral Hige;s bosons h? with masses rwA < ?Ti(h|') < me, where uia and 
me indicate the lower and upper bound of the mass spectrum. The squared coupling, g'^, of 
the Higgs states h? to the is modified by a factor fcj compared to the Standard Model H'^Z'^ 
couphng: g^o^o = h ■ ffzOhg^^ ■ 

If the Higgs states are assumed to be closer in mass than the experimental mass resolution, 
then there is no need to distinguish between separate ki. In this case the Higgs states and their 
reduction factors ki can be combined into a coupling density function, K{m) = dk/dm. The 
model obeys two sum rules which in the limit of unresolved mass peaks can be expressed as 
integrals over this coupling density function: 

oo 

J dm K{m) 



oo 

j dm K(m)m? 


where K{m) > and mc is a perturbative mass scale of the order of 200 GeV. The value of 
mc is model dependent and can be derived by requiring that there is no Landau pole up to a 
scale A where new physics occurs If neither a continuous nor a local excess is found in the 
data. Equation |2| can be used to place constraints on the coupling density function K{m). For 
example, if K{m) is assumed to be constant over the interval [m-A, m^] and zero elsewhere, 

K{m) = 1/ (mB — rnx) for ttt-a < m < me, 
= elsewhere, 

then certain choices for the interval [rriA, me] can be excluded. From this and from Equation ^ 
lower limits on the mass scale mc can be derived. 



(2) 



(3) 



2.2 The Stealthy Higgs scenario 

This scenario predicts the existence of additional SU(3)cxSU(2)lxU(1)y singlet fields (phions), 
which would not interact via the strong or electro-weak forces, thus coupling only to the Higgs 
boson Q. Therefore these singlets would reveal their existence only in the Higgs sector by offering 
invisible decay modes to the Higgs boson. The width of the Higgs resonance can become large 
if the number of such singlets, N, or the coupling u is large, thus yielding a broad spectrum in 
the mass recoiling against the reconstructed Z^. The interaction term between the Higgs and 
the additional phions in the Lagrangian is given by 

^interaction = "TTTtt'^^'/'V, (4) 

2v N 

where (p is the Standard Model Higgs doublet, u! is the coupling constant, and (p is the vector of 
the new phions. An analytic expression for the Higgs width can be found in the limit N ^ oo: 

FhC^Xh) = rsM("lH) + , (5) 

61 Vr TTlH 

where v is the vacuum expectation value of the Higgs field. This expression results when setting 
other model parameters to zero, including the mass of the phions The cross section for the 
Higgs-strahlung process can be calculated from Equations 9 and 10 of reference |^]. 
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In section ^.2.3 we derive limits on the Stealthy Higgs model which can be compared to 
expected limits from dedicated H*^ — > invisible searches, which are estimated in Ref. Q for the 
same scenario. By simulating signal spectra for different Higgs widths Fh we set limits in the 
uj-mii plane in the large N limit. 



3 Data sets and Monte Carlo samples 

The analyses are based on data collected with the OPAL detector at LEP during the runs in 
the years 1991 to 1995 at the peak (LEP 1) and on data taken in the years 1997 to 2000 at 
centre-of-mass energies between 183 and 209 GeV (LEP 2). The integrated luminosity used is 
115.4 pb~^ for the LEP 1 energy and 662.4 pb~^ for the LEP 2 energies, as detailed in Table ||. 
A description of the OPAL detector^ can be found elsewhere [||]. 

To estimate the detection efficiency for a signal from a new scalar boson and the amount of 
background from SM processes, several Monte Carlo samples are used. Signal events are simu- 
lated for masses from 1 keV to 110 GeV in a large variety of decay modes with the HZHA ||^] 
generator. The signal efficiencies are determined for all possible decays of a Standard Model Higgs 
boson (quarks, gluons, leptons, photons), for the decays into 'invisible' particles (e.g. Lightest 
Supersymmetric Particles) X^X^ as well as for 'nearly invisible' decays, S° X2X1) where 

the X2 decays into a Xi plus a photon or a virtual Z°, and for decays — > AA with A— > cc, 
gg or rr, where A is the CP-odd Higgs boson in supersymmetric extensions of the Standard 
Model. For simulation of background processes the following generators are used: BHWIDE ||^, 
TEEGG HI ((Z/7)* ^ e+e-(7)), KORALZ §, KK2F ||10| (both /U+Ai-(7) and r+r-(7)), JET- 
SET ||n|], PYTHIA Q (qq(7)), GRC4F ^ (four-fermion processes), PHOJET ||l|], HERWIG [|l|]. 



Vermaseren (hadronic and leptonic two-photon processes), NUNUGPV ||T6| (z^^'7) and RAD- 
COR ||l^ (77)- For all Monte Carlo generators other than HERWIG, the hadronisation is done 
using JETSET. The luminosity of the main background Monte Carlo samples is at least 4 times 
the statistics of the data for the two-fermion background, 50 times for the four-fermion back- 
ground and 5 times for the two-photon background. The signal Monte Carlo samples contain 
500-1000 events per mass and decay mode. The generated events are passed through a detailed 
simulation of the OPAL detector ||l8| and are reconstructed using the same algorithms as for the 
real data. 



4 Decay-mode independent searches for e+e ^S^Z^ 

The event selection is intended to be efficient for the complete spectrum of possible decay 
modes. As a consequence it is necessary to consider a large variety of background processes. 
Suppression of the background is performed using the smallest amount of information possible 
for a particular decay of the S^. The decays of the into electrons and muons are the channels 
with highest purity, and therefore these are used in this analysis. They are referred to as the 
electron and the muon channel, respectively. The signal process can be tagged by identifying 
events with an acoplanar, high momentum electron or muon pair. We use the term 'acoplanar' 
for lepton pairs if the two leptons and the beam axis are not consistent with lying in a single 
plane. 

Different kinematics of the processes in the LEP 1 and the LEP 2 analysis lead to different 
strategies for rejecting the background. At LEP 2 the invariant mass of the two final-state leptons 
in the signal channels is usually consistent with the Z^ mass, while this is not true for a large part 

^OPAL uses a right handed coordinate system. The z axis points along the direction of the electron beam and 
the X axis is horizontal pointing towards the centre of the LEP ring. The polar angle 9 is measured with respect 
to the 2 axis, the azimuthal angle <j) with respect to the x axis. 
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of the background. Therefore a cut on the invariant mass rejects a large amount of background. 
Remaining two-fermion background from radiative processes can partially be removed by using 
a photon veto without losing efficiency for photonic decays of the S^. In the LEP 1 analysis 
the invariant mass of the lepton pair cannot be constrained. Therefore, stronger selection cuts 
have to be applied to suppress the background, resulting in an insensitivity to the decays — > 
photons and at low masses also to — > e"'"e~. Hence, these decay modes are recovered in a 
search dedicated to e+e" — > S'^Z'^ with vv and S*^ — > photons (or photons plus invisible 

particles) or electrons at low masses mgo < 500 MeV. 

4.1 Event selection for e+e~ — > S°Z° with Z° — > e+e~ or ii^ ii~ 

The analysis starts with a preselection of events that contain at least two charged particles 
identified as electrons or muons. A particle is identified as an electron or muon, if it is identified 
by at least one of the two methods: 



The standard OPAL procedures for electron and muon identification |19|. These routines 
were developed to identify leptons in a hadronic environment. Since the signal events 
contain primarily isolated leptons, a second method with a higher efficiency is also used: 

A track is classified as an electron if the ratio Ej-p is greater than 0.8, where -p is the 
track momentum and E the associated electromagnetic energy. Furthermore the energy 
loss d£^/dx in the central tracking chamber has to be within the central range of values 
where 99 % of the electrons with this momentum are expected. Muons are required to have 
Ej-p < 0.2 and at least three hits in total in the muon chambers plus the last three layers 
of the hadronic calorimeter. 



The two tracks must have opposite charge and high momentum. Depending on the recoil mass 
of the lepton pair, the LEP 1 analysis requires a momentum of the higher energy lepton above 
20-27 GeV in the electron channel and above 20-30 GeV in the muon channel. The momentum 
of the lower energy lepton has to be greater than 10-20 GeV in both channels. 

For electrons these cuts apply to the energy deposited in the electromagnetic calorimeter, and 
for the muons to the momentum measured in the tracking system. At LEP 2 energies the lepton 
momenta have a weaker dependence on the recoil mass, therefore fixed cuts are used which are 
adjusted for the different centre-of-mass energies: Ei > 0.22 • y^, E2 > 0.11 • ^/s for electrons 
and pi > 0.22 • ^/s, p2 > 0.12 • ^/s for muons, where Ei, pi and E2, P2 are the energy and 
momentum of the lepton with the higher and lower momentum, respectively. 

The two leptons must be isolated from the rest of the event. The isolation angle aiso of a 
lepton candidate is defined as the maximum angle for which the energy -Econe contained within 
a cone of half-angle aiso around the direction of the lepton at the vertex is less than 1 GeV. 
^cone is the energy of all tracks and electromagnetic clusters not associated to a track within 
the cone, excluding the energy of the lepton itself. Leptons at small angles to the beam axis 
(I cos^l > 0.9 in the electron channel and | cos6\ > 0.94 in the muon channel) are not used due 
to detector inefficiencies and mismodelling in this region. These cuts also serve to reduce the 
background from two-fermion and two-photon processes. We ignore lepton candidates inside a 
0.3° azimuthal angle to the anode planes of the jet chamber since they are not well described in 
the detector simulation. If more than one electron or muon pair candidate with opposite charge 
is found, for the LEP 1 analysis the two leptons with the highest momentum, and for the LEP 2 
analysis the pair with invariant mass closest to mgo are taken as decay products. 

The background to the £~^£~ signal arises from several processes which are suppressed as 
described below: 
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• In (Z/7)* — > e"'"e~,/i+^^ events without initial or final state radiation the leptons are 
produced in a back-to-back topology. We reject these events by cutting on the acoplanarity 
angle (pa which is defined as vr — i;^>open, where (popen is the opening angle between the 
two lepton tracks in the plane perpendicular to the beam axis. For the LEP 1 analysis 
the acoplanarity angle is multiplied by the average of the sin 9 of the tracks in order to 
account for the larger influence of the track direction resolution on the acoplanarity angle 
for tracks with small sinO. The modified acoplanarity angle is termed a. The cuts are 
0.11 rad < a < 2.0 rad and (pa > 0.15 — 0.20 rad (depending on the centre-of-mass energy). 

• In two-photon processes, where the incoming electron and positron are scattered at low 
angles, usually one or both of the electrons are undetected. Events of this type usually 
have large missing momentum with the missing momentum vector, Pmiss, pointing at low 
angles to the beam axis. In (Z/7)* — > e~ , iJ,~^ fi~ events with initial-state radiation the 
photons usually remain undetected at low angles. The requirement | cos 0(pmiss)| < 0.98 
for the LEP 1 analysis and | cos 0(pmiss)| < 0.95 for the LEP 2 analysis reduces background 
from these two sources. 

• The semileptonic decays of b- or c-mesons provide another source of leptons which can be 
misidentified as direct decay products. This background is reduced by requiring the 
leptons to be isolated from the rest of the event. We require one of the isolation angles 
of the two lepton candidates to be greater than 20° and the other one to be greater than 
10° for the LEP 1 analysis and to be greater than 15° and 10°, respectively, for the LEP 2 
analysis. 

Up to this point the analyses for LEP 1 and LEP 2 energies are essentially identical, but they 
are tuned separately, as detailed in Table The different features of signal and background at 
LEP 1 and LEP 2 energies are taken into account with the following cuts. 



4.1.1 Cuts used only in the LEP 1 selection 

• Since the electron or muon pair originates from a its invariant mass is high in comparison 
to a typical pair of isolated leptons in hadronic background. We therefore require the lepton 
pair invariant mass to exceed 20 GeV. 

• At this stage the cut selection is still sensitive to all decay modes of the S". The main 
background, however, arises from electron and muon pairs accompanied by energetic photon 
radiation. Reduction of this kind of background is made by applying cuts on photons and 
electrons recognised as coming from a photon conversion. 

Events with less than four tracks are vetoed if there is an unassociated cluster in the 
electromagnetic calorimeter with an energy greater than 1 GeV outside a 10° cone around 
a lepton candidate (photon veto). They are also vetoed if the energy in the forward 
calorimeters, corresponding to the polar angle region 47-200 mrad, exceeds 2 GeV (forward 
veto). In order to reject events where the photon converts into an electron-positron pair, 
events with one, two or three tracks in addition to the lepton are excluded if at least one 
of them is identified as a track from a conversion (conversion veto). The conversion finder 
is based on an artificial neural network pO| . 

The photon and the conversion veto are at the expense of sensitivity for decays 
photons (or photons plus invisible particles) in the whole mass region and for — > e"'"e~ 
at low masses (m-go < 500 MeV). In order to retain sensitivity to these decay modes, they 
are taken into account in a search dedicated to e"'"e~ — > S'^Z" with Z^ — > uiy and —^ 
photons (or photons plus invisible particles) and for rrigo < 500 MeV also to electrons as 



described in section 4.2 
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All cuts are listed in Table § and the number of events after each cut is given in Table ^. The 
distributions of the cut variables in data and Monte Carlo are shown in Figures |and|. After the 
selection 45 events remain in the channel Z'^ — > e^e~, with 55.2 it 3.0 (stat.) it 3.0 (syst.) events 
expected from SM background (the evaluation of the systematic uncertainties is described in 
section 4.1.4). In the channel —>■ f^'^fJ.~ ^ 66 events remain in the data with 53.6 ± 2.7 (stat.) it 



2.1 (syst.) expected from SM background. 

Figure ^ shows the efficiency versus the mass for some example decay modes. The signal 
efficiency is at least 20% in the electron channel and at least 27% in the muon channel for 
masses between 4 and 45 GeV. At masses below the kinematic threshold for the decay of the 
into e"'"e~ (~ 1 MeV) only decays into photons or invisible particles are possible. For each 
mass hypothesis the smallest efficiency of all decay channels studied (also shown in Figure ^) is 
used in the limit calculation. The analysis is sensitive to a large range of masses, down to 
masses mgo well below T^o, where the cross section increases significantly. For this mass range 
mainly soft bosons with energy E^o < are emitted, but the spectrum exhibits a significant 
tail to large energies, which yields a detectable event topology. Figure ^ shows the recoil mass 
spectrum to the Z° decay products for both channels at ^/s = 91.2 GeV. The recoil mass squared 
is calculated from 



m 



rec 



{./s - Eu)'^ - p%, (6) 



where E^ and pu are the energy and the momentum sum of the two lepton tracks, and is 
the centre-of-mass energy. The momentum sum is calculated from the track momentum of the 
decay products in the muon channel and from the track momentum and energy deposition of 
the electrons in the electromagnetic calorimeter in the electron channel^. 

4.1.2 Cuts used only in the LEP 2 selection 

In the analysis for LEP 2 energies, signal and background characteristics differ significantly from 
those at LEP 1. 



• The most important difference compared to the LEP 1 analysis is the fact that in signal 
processes an on-shell boson is produced. The selection requires the invariant mass mu 
of the lepton pair to be consistent with the mass. Due to the limited detector mass 
resolution, invariant masses within mzo ± 8 GeV and mzo ± 10 GeV are accepted for the 
electron and the muon channel, respectively. 

• The dominant background at this stage originates from leptonic Z^ decays with photon 
radiation in the initial state. At centre-of-mass energies above m^o the cross section for 
radiating one (or more) high energy initial-state photon(s) is enhanced if the effective 
centre-of-mass energy of the electron-positron pair after photon emission is close to the Z" 
mass. Such events are called 'radiative returns' to the Z^ pole. These background events 
are characterised by an acolinear and sometimes acoplanar lepton pair and one or more 
high energy photons. Such events are rejected by a 7-veto: if there is only one cluster 
in the electromagnetic calorimeter not associated to a track and the energy -Eunass of the 
cluster exceeds 60 GeV, then the event is rejected. Events with two tracks and more than 
3 GeV energy deposition in the forward calorimeters (covering the polar angle region 47- 
200 mrad) are also vetoed. The cross section for two fermion production is much smaller 
at LEP 2 than at LEP 1 so events with final state radiation are not such an important 
background as in the LEP 1 case. 

^Due to the limited energy and momentum resolution, the calculated value of rriroc can be negative. We define 

rrircc ~ V'^roc foi m^ac ^ and mrcc = —\/—m'^cc for rn^ac < 0. 
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• In the remaining background from two-photon processes and (Z/7)* e"'"e^, with 
initial-state radiation the leptons carry considerable momentum along the beam axis. We 
reject these events by requiring |pf + < 50 GeV where pf are the z-components of the 
momentum of the two lepton candidates. 

All cuts are listed in Table |2| and the number of events after each cut is given in Table ^. 
The distributions of the cut variables in data and Monte Carlo are shown in Figures ^ and ^ for 
data taken at \/s= 183-209 GeV. A total of 54 events remain in the data of 183-209 GeV in the 
channel — > e"'"e~, with 46.9 ± 0.6 (stat.) ± 3.3(syst.) events expected from SM background 
(the evaluation of the systematic uncertainties is described in section [1.1.41 ). In the channel 

—>■ fi'^fi", 43 events remain in the data with 51.6 ± 0.3 (stat.) ± 2.6 (syst.) expected from SM 
background. The signal efficiency is at least 24% in the electron channel and at least 30% in the 
muon channel for masses between 30 and 90 GeV. 

Figure |^ shows the efficiency versus the S° mass at ^/s = 202-209 GeV for some example 
decays as well as the minimum efficiencies which are used in the limit calculation. The efficiencies 
for 183-202 GeV have similar values for mgo < 100 GeV. For the lower centre-of-mass energies 
the efficiency decreases faster for higher masses due to kinematic effects, primarily the cut on 
the acoplanarity angle. Figure § shows the recoil mass spectrum for both channels summed from 
183-209 GeV. 

4.1.3 Correction on background and signal efficiencies 

In all channels a correction is applied to the number of expected background events and the 
signal efficiencies due to noise in the detectors in the forward region which is not modelled by 
the Monte Carlo. The correction factor is derived from the study of random beam crossings. 
The fraction of events that fail the veto on activity in the forward region is 7.5 % for LEP 1 and 
3.1 % for LEP 2. Since the veto is applied only to events with less than or equal to four tracks, 
the corrections on the expected background in the actual analyses are typically only 1.8-3.5%. 
For the signal efficiencies the full correction is applied to the decay channels where appropriate. 

4.1.4 Systematic uncertainties 

The systematic uncertainty of the lepton identification efficiency is studied in a control sample 
of events with two collinear tracks of which at least one is tagged as an electron or muon. The 
systematic uncertainty is obtained from the difference of the identification efficiencies for the 
other track between data and Monte Carlo. 

The tracking systematics are studied by changing the track resolution^ in the Monte Carlo 
by a relative fraction of 5% in do and cj) and by 10% in zq, k and cot 6, which corresponds to the 
typical difference in the resolution of these parameters in data and Monte Carlo. The difference 
in signal and background expectation compared to the one obtained from the unchanged track 
resolution is taken as the systematic uncertainty. 

The reconstruction of the energy deposition in the electromagnetic calorimeter and the mo- 
mentum in the tracking system of the lepton candidates is investigated with the help of the mean 
values Xdat and xmc of the distributions of p and E from the collinear lepton pair control sample 
for data and Monte Carlo expectation. The analyses are repeated with the cuts on p and E 
being changed by the difference \xdat — xmc\- The deviations in the number of expected events 
compared to the original cuts are taken as the systematic uncertainties. 

The uncertainty from the lepton isolation angle Oiso is studied in different ways for the LEP 1 
and LEP 2 analysis. In the LEP 1 selection a control sample of hadronic events is selected. 

*do is the distance between the vertex and the point of closest approach of a track to the vertex in the r-(l) 
plane, 20 is the z-coordinate of the track at this point, and k is its curvature. 
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Random directions are then chosen in the event, and the angles aiso of the vectors pointing 
to these directions are determined. The cut on aiso is varied by the difference of the mean of 
the data and the Monte Carlo distributions of the control sample. In the LEP 2 selection the 
uncertainty is obtained in a similar way but from the isolation of the lepton in W'^'W" qqli' 
events. 

Correct modelling of photon radiation and conversions is a crucial ingredient of the decay- 
mode independent searches. For LEP 1, the effect of the description of photon radiation in 
the Monte Carlo is estimated from the difference in the number of events between data and 
background expectation after removing the photon and conversion veto. At least one identified 
photon or conversion is required for the tested events. In the muon channel at LEP 2 energies 
two different Monte Carlo generators are used for the two-fermion background, and the difference 
between the background prediction of the two generators is taken as the systematic uncertainty 
of the photon modelling. For the electron channel only one generator is available. Here, the 
uncertainty is determined from the comparison of the number of events in the data and Monte 
Carlo sets in a side band of the distribution of the lepton pair invariant mass where no signal is 
expected. This test is dominated by the statistical uncertainties of the side-band sub-sample. 

In the analysis the four-fermion Monte Carlo samples are reweighted to account for low mass 
resonances (e.g. p, uj, (p, J ftp) and the running of aem(9^)- The uncertainty from this reweighting 
is assessed to be 50% of the change of the expected background after switching off the reweighting. 

All uncertainties for a particular centre-of-mass energy are assumed to be uncorrelated and 
the individual contributions are added in quadrature for the total systematic uncertainty. The 
dominant systematic uncertainties in the LEP 1 background expectation come from the de- 
scription of photon radiation and photon conversions in the Monte Carlo as well as from the 
uncertainty of the four-fermion cross section. The precision of the predicted signal efficiency is 
mainly limited by the description of the lepton isolation. 

In the LEP 2 selection the modelling of the radiative returns has a large impact on the total 
systematic uncertainty, both in the electron and the muon channel. In the electron channel the 
uncertainty from the isolation angle criterion and in the muon channel the uncertainty of the 
muon identification efficiency are also significant. 

For the LEP 2 data, the evaluation of the systematic uncertainties at each single centre-of- 
mass energy is limited by Monte Carlo statistics. Therefore they are investigated for the total 
set of Monte Carlo samples with = 183-209 GeV. 

The numbers of expected background events for the LEP 1 and LEP 2 analyses, broken down 
by the different centre-of-mass energies, are listed in Table ^ for the channels — > e"'"e~ and 
fj-^fj,^. The numbers include systematic errors discussed above and the statistical error from 
the limited Monte Carlo samples. Also the number of expected events from a 30 GeV Standard 
Model Higgs boson is shown. A detailed overview of the different systematic uncertainties is 
given m Table |. 

4.2 Event selection for e+e~ — ^ S°Z° — ^ n • 7 uT', e^e~vV 

In this section a search for e"'"e~ — > S'^Z'^ n-'^vV ox e^e~i>V (the latter for mgo < 500 MeV) at 
y/s = 91.2 GeV is described. The removal of radiative backgrounds in the selection described in 
section [4.1.1| rejects the S*^ decays into photons (due to the photon veto) and, in the mass region 
mgo < 500 MeV, also the decays into electrons (due to the conversion veto). The selection for 
S'^Z'^ S^vV is included to recuperate the sensitivity to the photon and electron decay modes 
and therefore to remain decay-mode independent. 
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4.2.1 Event selection 



In signal events the S° is radiated off the TP with some kinetic energy and a certain amount of 
transverse momentum. Therefore the total visible energy i?vis in the electromagnetic calorimeter 
is required to exceed 12 GeV and the transverse momentum pT reconstructed from the event 
is required to exceed 6 GeV. Since the decays into neutrinos which carry energy out of the 
detector, the total amount of visible energy is reduced. The selection requires -Evis < 60 GeV. 

Several cuts are applied to reduce the background from processes with topologies different 
from the signal. The selection allows only events with zero {n ■ 7 zvI7-channel) or two (e"'"e~ 
I'F-channel) identified electrons, using the same electron identification routines as described in 
section [l.l| . The next selection cuts are intended to further reduce background from cosmic rays 
or beam halo particles. Events triggered by cosmic rays or beam halo particles are characterised 
by extended clusters in the electromagnetic calorimeter, hits in the hadron calorimeter and muon 
chambers and a signal from the time-of-flight counter that shows a significant discrepancy from 
its expected value. We therefore require no hit in the muon chambers and at most two hits in 
the hadron calorimeter. No more than one cluster with an energy deposition larger than 2 GeV 
is allowed in the hadron calorimeter. The number of lead glass blocks in each cluster in the 
electromagnetic calorimeter must be less than 15. The difference between the measured time of 
flight and the expected time for a particle coming from the interaction point is required to be 
less than 2 ns. 

The remaining background is mostly from e"'"e~ — > events, where the photon is usually 
emitted at small angles to the beam axis. A hard cut on the angular distribution of clusters in 
the electromagnetic calorimeter and the forward calorimeters (polar angle region 47-200 mrad) 
is applied. For this purpose the polar angle of the energy vector E is defined as 



E 



Ei ■ 9i 



■ , E^is 

The sum runs over all clusters i (with polar angle 9i). Cuts applied on the energy vector and the 
individual clusters are | cos0^| < 0.65 and | 0036*41 < 0.9. The energy in the forward calorimeters 
must be less than 2 GeV. 

Events with two electrons must satisfy some additional requirements. The tracks must be 
identified as electrons with opposite charge. The angle A(/) between the tracks must be less than 
10°, the invariant mass niu must be less than 2 GeV, and a transverse momentum pT > 7 GeV 
of the event is required. Events are rejected if there are any additional clusters other than those 
associated with one of the two tracks. 



A correction due to random detector occupancy is applied as described in section ^.1.3| . The 
full correction of —7.5 % is used since the forward detector veto applies to all events. 

After all cuts 15 events are selected from the data with a background prediction of 11.3 it 
l.l(stat.) lb 0.2(syst.), where the uncertainties are evaluated as described below. Figure^ shows 
the distribution of the cut variables in data and Monte Carlo. 

There is no statistically significant excess in this channel, and the shape of the distributions 
of the cut variables in data and Monte Carlo are in good agreement. Furthermore the product of 
the signal efficiency and the — > vu decay branching ratio is substantially higher than for the 

^'^^"-channels, and the predicted background is much less. Hence, the channel — > uu has 
much higher sensitivity than the electron and the muon channel. It does not contribute to the 
actual limits, provided that the systematic uncertainties are not much larger than in the other 
channels. For a conservative limit calculation only the channels with lowest sensitivity are used. 

The search channel described in this section recovers sensitivity to the decay modes n-j 
and at low masses mgo to the decay S" — > e+e" to which the analysis described in section |4.1.1] 
had no sensitivity. However, the requirement E'vis > 12 GeV in this channel can lead to an 
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insensitivity to decays — > X2X1 ~^ 7XiXi for certain combinations of mgo,m^o and m^o for 
the whole LEP 1 analysis. 

4.2.2 Systematic uncertainties 

Uncertainties in this channel predominantly come from the energy calibration of the electromag- 
netic calorimeters. In reference |^ it is shown that an electromagnetic cluster has a calibration 
uncertainty of 25 MeV. Since the number of clusters in the electromagnetic calorimeter for se- 
lected data and Monte Carlo events is less than five, the deviation of the expected number of 
background events and the signal efhciencies after shifting the cuts on the visible energy i?vis 
and the transverse momentum pT by four times 25 MeV is used as a systematic uncertainty. 



The deviation is found to be 1.4%. Prom the same reference |^ we take the systematic un- 
certainties on the time-of-flight signal (0.5%). Other sources for systematic uncertainties are 
the luminosity (0.5%), the limited Monte Carlo statistics (10.0% for background and 2.2% for 
the signal) and, for events where the decays into electrons, the uncertainty on the electron 
identification efficiency (0.8%). Summing these individual sources up in quadrature, estimates 
of the total uncertainty in the background of 10.0 %(stat.)-|-1.8 %(syst.) and in the signal of 
2.2 %(stat.)-|-1.6 %(syst.) are obtained. Given the expected number of signal and background 
events, this is much less than the level where the channel starts contributing to the limit calcu- 
lation. 



5 Results 

The results of the decay-mode independent searches are summarised in Table ||, which compares 
the numbers of observed candidates with the background expectations. The total number of 
observed candidates from all channels combined is 208, while the Standard Model background 
expectation amounts to 207.3±4.1(stat.) it ll.l(syst.). For each individual search channel there 
is good agreement between the expected background events and observed candidates. As no 
significant excess over the expected background is observed in the data, limits on the cross 
section for the Bjorken process e"'"e~ S^Z^ are calculated. 

The limits are presented in terms of a scale factor k, which relates the cross section for S'^Z'^ to 
the Standard Model one for the Higgs-strahlung process e+e" — > Hgj^Z'' as defined in Equation ||. 
The 95% CL upper bound on k is obtained from a test statistic for the signal hypothesis, by 



using the weighted event-counting method described in |^: In each search channel, given by the 
different centre-of-mass energies and the Z" decay modes considered, the observed recoil mass 
spectrum is compared to the spectra of the background and the signal. The latter is normalised 
to £ ■ BR ■ C ■ k ■ cthozo , where e is the minimum signal detection efficiency out of all tested 
decay modes, BR is the branching ratio of the decay mode considered in this channel and C 
is the integrated luminosity recorded for that channel. The efficiencies for arbitrary masses 
are interpolated from the efficiencies at masses for which Monte Carlo samples were generated. 
Every event in each of these mass spectra and each search channel is given a weight depending 
on its expected ratio of signal over background, s/b, at the given recoil mass. For every assumed 
signal mass these weights are a function of the signal cross section, which is taken to be k 
times the Standard Model Higgs cross section for the same mass. Finally, from the sum of 
weights for the observed number of events, an upper limit k^^ for the scale factor is determined 
at the 95% confidence level. 

The systematic uncertainties on the background expectations and signal selection efficiencies 



are incorporated using the method described in |23| 
The limits are given for three different scenarios 

1. Production of a single new scalar S°. 
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2. The Uniform Higgs scenario. 



3. The Stealthy Higgs scenario. 
5.1 Production of a single new scalar S° 

In the most general interpretation of our results, a cross-section limit is set on the production of 
a new neutral scalar boson S° in association with a boson. To calculate the limit we use the 
mass distributions of which the sums are shown in Figure ^ and |8| for OPAL data, the expected 
background and the signal. 



In Figure |11| we present the limits obtained for scalar masses down to the lowest generated 
signal mass of 1 keV. They are valid for the decays of the S'' into hadrons, leptons, photons and 
invisible particles (which may decay inside the detector) as well as for the case in which the S" 
has a sufficiently long lifetime to escape the detector without interacting or decaying. A decay of 
the S'^ into invisible particles plus photons, however, can lead to a reduced sensitivity in the mass 



region where the sensitivity of the analyses is dominated by the LEP 1 data (see section 4.2. 1| ). 



The observed limits are given by the solid line, while the expected sensitivity, determined from 
a large number of Monte Carlo experiments with only background, is indicated by the dotted 
line. The shaded bands indicate the one and two sigma deviations from the expected sensitivity. 
Values of > 0.1 are excluded for values of mgo below 19 GeV, whereas /c > 1 is excluded from 
the data for mgo up to 81 GeV, independent of the decay modes of the S° boson. This means 
that the existence of a Higgs boson produced at the SM rate can be excluded up to this mass 
even from decay-mode independent searches. For masses of the new scalar particle well below 
the width of the Z°, i.e. mgo < 1 GeV, the obtained limits remain constant at the level of 
A;95^ = 0.067, and /c^Sp = 0.051. 

The discrepancy between the expected and the observed limits is within one standard devia- 
tion for masses below 52 GeV and for masses above 82 GeV. The deviation of about two sigma in 
the mass range 52-82 GeV is due to a deficit of selected data events in the recoil mass spectrum 
of both the electron and muon channels. 

5.2 Limits on signal mass continua 
5.2.1 The Uniform Higgs scenario 

We simulated signal spectra for the Uniform Higgs scenario for K = constant over the interval 
[ttt-a, ?7t,b] and zero elsewhere. Both the lower mass bound mA and the upper bound me are 
varied between 1 GeV and 350 GeV (with the constraint ttia < me). In a similar way to the 
previous section we get an upper limit on the integral in Equation |2|. 

Figure shows the mass points (mp^^m^) for which the obtained 95% CL limit on J dmK 
is less than one. These are the signal mass ranges m\ < mjjp < me which can be excluded 

assuming a constant K. 

If ruA = rriB, then the signal spectrum reduces to the mass distribution of a single boson. 



Excluded points on the diagonal niA = tub are therefore the same masses as in Figure 11 for 
which k < 1. The horizontal line illustrates an example for excluded mass ranges: The line starts 
on the diagonal at mA = tub = 35 GeV and ends at uib = 99 GeV. This value of tub is the 
highest upper mass bound which can be excluded for this value of tua- All mass ranges with an 
upper bound me below 99 GeV are also excluded for ttia = 35 GeV. The highest excluded value 
of ruB {rriB = 301 GeV) is achieved for uia set to GeV. 



Using the two sum rules from section lower limits on the perturbative mass scale mc 



can be derived. For each excluded value of uia we take the highest excluded value of uib and 
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determine the lower bound of mc according to Equation y. The excluded mass ranges for mc, 
assuming a constant K, are shown in Figure [l^. 

5.2.2 Bin-by-bin limits 



The limits presented in section |5.2.1| are specific to the case where the coupling density is constant 
in the interval [mA, ttib] and zero elsewhere. The data can also be used to exclude other forms of 
K{m). To provide practical information for such tests, we have measured K{m) in mass bins with 
a width comparable to the experimental mass resolution. The typical resolution of the recoil mass 
in the LEP 1 analysis varies between 1 and 5 GeV in the mass region between 10 and 55 GeV. In 
the LEP 2 analysis the width is between 3 and 15 GeV for recoil masses between 20 and 100 GeV. 
The width gets smaller at higher recoil masses. The results of the measurement of K{m) are 
given in Table ^ together with the corresponding statistical and systematic uncertainties. 

From these measured numbers oiK{m)^ one can obtain upper limits on the integral J dmK{m) 
for any assumed shape of K{m) using a simple fitting procedure. To account for mass resolu- 
tion effects, we provide a correction matrix (Table P). To test a certain theory with a distribution 

of K{m) values in the 10 measured bins from Table |^, written as a vector k = (^Ki, . . . , Kio^ , 

the corrected vector k' = Ck, can then be fitted to the measured values. In the fit the systematic 
uncertainties, which are small compared to the statistical errors, can be assumed to be fully 
correlated bin-by-bin. 

5.2.3 The Stealthy Higgs scenario 

To set limits on the Stealthy Higgs scenario we have simulated the spectrum of a Higgs boson 
with a width according to Equation |5| and Ref . Q . 



The excluded regions in the uj-mn parameter space are shown in Figure |14[ To illustrate 
the Higgs width according to Equation ^, for a given mass ttt-h and coupling u> 'isolines' for 
some sample widths are added to the plot. The vertical edge in the exclusion contour at mn = 
81 (62) GeV in the observed (expected) limits refiects the detector mass resolution in Am: For a 
fixed mass mn the exclusion power is the same for all couplings u; that yield Fh < Am, and the 
limits for u; ^ reproduce the limits for a single narrow in Figure |ll[ The maximal excluded 
region of the coupling u is achieved for masses around 30 GeV, where to can be excluded up to 
Lj = 2.7. For lower masses the sensitivity drops due to the rapidly increasing width of the Higgs 
boson, and for higher masses due to the decreasing signal cross section. 



6 Conclusions 

Searches for new neutral scalar bosons S° decaying to hadrons of any fiavour, to leptons, photons 
invisible particles and other modes have been performed based on the data collected at ^/s = 
mz and 183 to 209 GeV by studying the recoil mass spectrum of e+e", ;U+/u~ in S'^Z'^ 

production and the channel where the decays into uu and the into photons or e"'"e~. No 
significant excess of candidates in the data over the expected Standard Model background has 
been observed. Therefore upper limits on the production cross section for associated production 
of and Z", with arbitrary decay modes, were set at the 95 % confidence level. Upper limits in 
units of the Standard Model Higgs-strahlung cross section of /c < 0.1 for 1 keV < mgo < 19 GeV 
and /c < 1 for mgo < 81 GeV were obtained. In further interpretations, limits on broad continuous 
signal mass shapes to which previous analyses at LEP had no or only little sensitivity were set for 
the first time. Two general scenarios in the Higgs sector were investigated: A uniform scenario, 
when the signal arises from many unresolved Higgs bosons, and a Stealthy Higgs model, when 
the Higgs resonance width is large due to large Higgs-phion couplings. 
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\/s (GeV) 


year 


integrated luminosity (pb ^) 


91.2 


1989-95 


115.4 


183 


1997 


56.1 


189 


1998 


177.7 


192 


1999 


28.8 


196 


1999 


73.2 


200 


1999 


74.2 


202 


1999 


36.5 


202-206 


2000 


83.1 


206-209 


2000 


132.4 



Table 1: Overview of the analysed integrated data luminosities. 



LEP 1: ^ e+e-,At+/i- 


0. Preselection 


see text 


1. Modified acoplanarity 


0.11 rad < a < 2.0 rad 


2. Polar angle of missing momentum vector 


1 cos6i(pmiss)| < 0.98 for Pmi^s > 2 GeV 


3. Isolation of lepton tracks 


max(aisoi,aiso2) > 20° 




min(aisoi,aiso2) > 10° 


4. Invariant mass of the lepton pair 


20 GeV < mu < 100 GeV 


5. Photon and Conversion veto 


see text 



LEP 2: ^ e+e-,Ai+M" 


0. Preselection 


see text 


1. Acoplanarity 


(t>a > 0.15-0.20 rad 


2. Polar angle of missing momentum vector 


1 cos6'(pmiss)| < 0.95 for pmiss > 5 GeV 


3. Isolation of lepton tracks 


max(Q;isoi,aiso2) > 15° 




min(aisoi,aiso2) > 10° 


4. Invariant mass of the lepton pair 


TTlg+e- — TOzo < 8 GeV 




- "^zo| < 10 GeV 


5. Photon veto 


see text 


6. Momentum in ^-direction 


\pi < 50 GeV 



LEP 1: Z" ^ i/z7 


1. Cosmic muon and beam halo veto 


see text 


2. Number of identified electron tracks 


A^e = or 2 


3. Visible energy in electromagnetic calorimeter 


E^Ecai > 12 GeV, < 60 GeV 


4. Transverse momentum of event 


PT > 6 GeV 


5. Direction of energy vector 


1 cos 61^1 < 0.65 


6. Energy in forward detector 


EFdet < 2 GeV 


Additional cuts for events with two electron tracks 




7. Angle between tracks 


< 10° 


8. Transverse momentum of event 


Pt > 7 GeV 


9. Unassociated clusters in electromagnetic calorimeter 


N — 

-^'unass. — ^ 



Table 2: A summary of the selection criteria. 
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= 91.2 GeV 








Cut 


Data 


Total 


2-fermion 


4-fermion 


2-photon 


Signal 






bkg. 






(mso=30 GeV) 


Electron channel 


Preselection 


122431 


129115 


128490 


586.3 


38.9 


46.2 % 


a 


1560 


1694 


1628 


58 


8 


33.4% 




1500 


1628 


1571 


55 


2 


32.8% 


Lepton isolation 


1368 


1466 


1414 


50 


2 


28.6% 


Mu 


1362 


1462 


1410 


50 


2 


28.6% 


Photon+ Conversion veto 


45 


55.2 


20.5 


34.4 


0.3 


28.6% 


Muon channel 


Preselection 


109552 


115001 


114475 


459.1 


66.6 


54.0 % 


a 


1575 


1601 


1526 


58 


17 


40.2 % 


|cos6ip„.^J 


1549 


1575 


1512 


57 


6 


40.0% 


Lepton isolation 


1403 


1470 


1412 


52 


6 


37.4% 


Mu 


1397 


1467 


1410 


51 


6 


37.4% 


Photon+ Conversion veto 


66 


53.6 


17.0 


35.4 


1.2 


35.0% 



Cut 


Data 


Total 
bkg. 


ui"y 


leptons 


other 


Signal 


Missing energy channel 




Events with tracks 






mso=5 GeV, — > 77 


Preselection 




73 


68.5 


63.5 


4.8 


0.3 


44.2 % 


1 cos 6*^1 < 0.65 




54 


51.1 


48.1 


2.7 


0.3 


38.4% 


^Ecai > 12 GeV 




14 


10.7 


9.8 


0.6 


0.3 


30.0% 




Events with 2 tracks 






mso=100 MeV, ^ e+e" 


Preselection 




30 


21.6 


4.5 


11.6 


5.5 


29.2% 


1 cos 6*^1 < 0.65 




17 


14.4 


3.7 


9.5 


1.2 


25.7% 


Acf) < 10° 




13 


7.9 


3.6 


3.7 


0.6 


25.7% 


rn-e+c- < 2 GeV 




12 


7.9 


3.6 


3.7 


0.6 


25.7% 


Charge qi ■ q2 = - 


1 


10 


6.3 


3.6 


2.1 


0.6 


25.2% 


■^unass. — 




4 


4.0 


3.4 


0.0 


0.6 


23.7% 


> 7 GeV 




3 


2.5 


2.5 


0.0 


0.0 


20.9% 


^Ecai > 12 GeV 




1 


0.5 


0.5 


0.0 


0.0 


14.8 % 



Table 3: Cutflow tables for the LEP 1 analysis: Number of selected events after each cut. As an 
example the efficiencies for the signal process with S° —>■ hh are given for the lepton channels, and 
with S° — > 77 and — > e"'"e~ for the missing energy channel. The preselection in the missing 
energy channel includes the cuts 1, 4, 6 from the vV channel in Table ^, -Etotai > 4 GeV 

and I cos^il < 0.9 (see Section 4.2.1 ). 
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= 183-209 GeV 


Cut 


Data 


Total 


2-fermion 


4-fermion 


2-photon 


Signal 






bkg. 






(toso=90 GeV) 


Electron channel 


Preselection 


27708 


28183.5 


27720.0 


378.0 


85.0 


49.1% 


Lepton isolation 


24176 


24803.9 


24410.6 


314.3 


79.0 


42.1% 


Mee 


708 


639.1 


547.9 


73.0 


18.2 


37.7% 


Photon-veto 


470 


477.1 


393.8 


67.9 


15.4 


37.7% 


1 cos 6'p„i^ J 


118 


106.3 


57.4 


45.7 


3.2 


34.8 % 


Acoplanarity 


67 


63.1 


25.4 


37.2 


0.5 


28.7% 




54 


46.9 


12.8 


33.7 


0.4 


28.7% 


Muon channel 


Preselection 


3042 


3115.6 


2818.8 


212.2 


84.6 


64.7% 


Lepton isolation 


2866 


2948.5 


2669.5 


195.9 


83.1 


55.7% 


Mee 


803 


842.4 


733.3 


88.5 


20.7 


49.3 % 


Photon-veto 


575 


629.3 


532.0 


80.9 


16.4 


49.3 % 


|cos6'p,„,J 


111 


101.5 


45.8 


52.3 


3.4 


45.5 % 


Acoplanarity 


66 


72.0 


26.7 


44.3 


1.0 


37.5 % 


bf +P2I 


43 


51.6 


12.2 


38.6 


0.8 


37.5 % 



Table 4: Cutflow tables for the LEP 2 analysis: Number of selected events after each cut. As 
an example the efHciencies for the signal process S°Z° — bb^''"^" are also given. The efficiencies 
are the average of the values at 183-209 GeV. 
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(GeV) 


Data 


Total 


2-fermion 


4-fermion 


2-photon 


Signal 




bkg. 








(mso=30 GeV) 


Electron channel 


91.2 


45 


55.2±3.0±3.0 


20.5 


34.4 


0.3 


15.61 ±0.31 ±0.47 


183 


7 


3.6 ±0.1 ±0.3 


1.4 


2.1 


0.1 


0.91 ±0.02 ± 0.03 


189 


18 


13.7±0.4± 1.0 


4.2 


9.5 


0.0 


2.42 ±0.04 ±0.09 


192 





2.2 ±0.1 ±0.2 


0.7 


1.5 


0.0 


0.37±0.01±0.01 


196 


6 


5.7±0.2±0.4 


2.0 


3.7 


0.0 


0.87±0.01±0.03 


200 


4 


4.8 ±0.2 ±0.3 


1.2 


3.5 


0.1 


0.81 ±0.01 ±0.03 


202 


5 


2.5 ±0.1 ±0.2 


0.6 


1.9 


0.0 


0.39 ±0.01 ±0.01 


202-206 


5 


5.0 ±0.2 ±0.4 


0.7 


4.2 


0.1 


0.86 ±0.01 ±0.03 


206-209 


9 


9.4 ±0.3 ±0.7 


2.0 


7.3 


0.1 


1.34 ±0.02 ±0.05 


E(> 183) 


54 


4C.9±().(i± 3.5 


12.8 


33.7 


0.4 


7.97 ±().()(i± 0.25 


Muon channel 


91.2 


66 


53.6±2.7± 2.1 


17.0 


35.4 


1.2 


21.55 ±0.45 ±0.69 


183 


5 


4.4±0.1±0.2 


1.6 


2.7 


0.1 


1.20 ±0.01 ±0.05 


189 


9 


13.7±0.1±0.7 


4.0 


9.5 


0.2 


2.96 ±0.03 ±0.11 


192 


2 


2.5 ±0.1 ±0.1 


0.6 


1.9 


0.0 


0.46 ±0.01 ±0.02 


196 


6 


6.1 ±0.1 ±0.3 


1.2 


4.7 


0.2 


0.96 ±0.01 ±0.04 


200 


5 


5.7±0.1±0.3 


1.3 


4.3 


0.1 


0.89 ±0.01 ±0.03 


202 


3 


2.9 ±0.1 ±0.1 


0.6 


2.3 


0.0 


0.43 ±0.01 ±0.02 


202-206 


9 


6.0 ±0.1 ±0.3 


0.9 


5.0 


0.1 


1.00±0.01±0.04 


206-209 


4 


10.3 ±0.1 ±0.5 


2.0 


8.2 


0.1 


1.53 ±0.02 ±0.06 


E(> 183) 


43 


51.6±0.3±2.5 


12.2 


38.6 


0.8 


9.43 ±0.06 ±0.37 



Missing energy channel 


(GeV) 


Data 


Total 
bkg. 




n ■ 7 


other 


Signal 
(mso=30 GeV) 


91.2 


15 


11.3±1.1±0.2 


10.3 


0.3 


0.7 


175.07 ±3.85 ±2.80 



Table 5: Summary of selected data events, background Monte Carlo and signal expectation for 
a 30 GeV Standard Model Higgs boson in the decay-mode independent searches. The first error 
is statistical and the second error is systematic. 
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Electron channel - 


uncertainties in % 




Source 


91 GeV 


183-209 GeV 


Bkg. 


Sig. 


Bkg. 


Sig. 


Electron-ID 


0.8 


0.8 


1.3 


1.3 


Energy 






1.2 


1.5 


Isolation angle 


1.3 


2.6 


4.3 


2.8 


Trk. resolution 


2.3 


1.3 


2.2 


1.3 


ISR/FSR 


2.4 




4.7 






4.0 




0.4 




Luminosity 


0.5 




0.2 




Total systematics 


5.4 


3.0 


7.0 


3.7 


Statistics 


5.5 


2.0 


3.1 


0.9 



Muon channel - uncertainties in % 


Source 


91 GeV 


183-209 GeV 


Bkg. 


Sig. 


Bkg. 


Sig. 


Muon-ID 


1.5 


1.5 


2.8 


2.8 


Momentum 






1.9 


1.3 


Isolation angle 


0.2 


2.1 


1.7 


2.0 


Trk. resolution 


2.7 


1.9 


2.2 


1.1 


ISR/FSR 


2.0 




2.3 






1.2 




0.2 




Luminosity 


0.5 




0.2 




Total systematics 


3.9 


3.2 


5.0 


3.8 


Statistics 


5.1 


2.1 


1.7 


1.0 



Table 6: Systematic uncertainties in percent for background and signal. For ^/s = 91 GeV the 
uncertainties are given for mgo = 30 GeV, for ^/s > 91 GeV they are shown for mgo = 60 GeV. 





Measurement of K in bins of 10 GeV width 


Bin 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Mass (GeV) 


0-10 


10-20 


20-30 


30-40 


40-50 


50-60 


60-70 


70-80 


80-90 


90-100 


K 


2.1 


-2.4 


-4.9 


-2.8 


-7.1 


5.8 


-33.5 


-45.2 


-18.6 


200.2 


xl03 (GeV"^) 






















A(i^),tat. 


2.9 


4.4 


4.8 


6.4 


14.3 


23.7 


21.7 


30.5 


66.9 


166.4 


xl03 (GeV^^) 






















A(^)sys. 


0.9 


0.9 


0.6 


0.5 


1.3 


3.3 


4.5 


7.0 


16.5 


37.4 


xl03 (GeV-^) 























Table 7: Bin-wise measurement of K for the mass range 0-100 GeV with Am = 10 GeV. To 
fit a theoretical distribution K to these values, the correction matrix C from Table |^ must be 
applied first. 
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c = 



0.33 





04 





02 





01 





00 





00 





00 





00 


0.00 





00 


0.41 





53 





02 





01 





00 





00 





00 





00 


0.00 





00 


0.17 





27 





11 





04 





01 





00 





00 





00 


0.00 





00 


0.09 





11 





14 





17 





05 





01 





00 





00 


0.00 





00 


0.00 





05 





17 





21 





28 





09 





01 





00 


0.00 





00 


0.00 





00 





16 





20 





27 





34 





11 





01 


0.00 





00 


0.00 





00 





11 





14 





16 





29 





37 





09 


0.01 





00 


0.00 





00 





10 





10 





10 





15 





25 





43 


0.06 





00 


0.00 





00 





09 





07 





07 





08 





11 





27 


0.46 





03 


0.00 





00 





09 





07 





05 





04 





06 





11 


0.32 





32 



Table 8: Correction matrix for mass resolution. For a given theory to be tested with a distri- 
bution of K values in the 10 mass bins, k = (^Ki, . . . ,Kio^ , the vector k has to be multiplied 

by the matrix C to account for mass resolution effects. The corrected vector k' = Ck can then 
be fitted to the measured values of K from Table 0. 
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maxia.^^^,a.^J (deg) 
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20 
o 
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I I. ■ ■ L ■ ■ L' 1' L'^jf - i" 1 




50 100 150 

i- OPAL data 
I 4 fermion 
I 2 fermion 



mg = 1.2 GeV x 0.1 
mg = 30 GeV x 1.5 



20 40 60 80 100 

m,, (GeV) 

Figure 1: Cut variables for — > e"'"e~ at ^/s = 91.2 GeV. The OPAL data are indicated by dots 
with error bars (statistical error), the four-fermion background by the light grey histograms and 
the two-fermion background by the medium grey histograms. The signal distributions from a 
1.2 GeV are plotted as dashed lines and those from a 30 GeV as dotted lines, respectively. 
The signal histograms are normalised corresponding to 0.1 and 1.5 times of the Standard Model 
Higgs-strahlung cross section and show the decay channel — > gg- Each variable is shown with 
the cuts applied before the cut on this variable is done, respecting the order of cuts in Table ||. 
The arrows indicate the accepted regions. The histograms in Figure a) have non-constant bin 
widths. 
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Figure 2: Cut variables for A*^/^" at ^/s = 91.2 GeV. The OPAL data are indicated by 

dots with error bars (statistical error), the four- fermion background by the light grey histograms 
and the two-fermion background by the medium grey histograms. The signal distributions from 
a 1.2 GeV are plotted as dashed lines and those from a 30 GeV S" as dotted lines, respectively. 
The signal histograms are normalised corresponding to 0.1 and 1.5 times of the Standard Model 
Higgs-strahlung cross section and show the decay channel — > gg- Each variable is shown with 
the cuts applied before the cut on this variable is done, respecting the order of cuts in Table ||. 
The arrows indicate the accepted regions. The histograms in Figure a) have non-constant bin 
widths. 
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Figure 3: The efficiency versus the mass at -y/s = 91.2 GeV for a subset of decay modes 
of S^: a)+b) — > e"'"e~ in hnear and logarithmic mass scale; c)+d) —>■ fi^ iJi~ in linear and 
logarithmic mass scale. The minimum efficiencies which are used in the limits are also shown. In 
the low mass region, below the threshold for the decays of the S*^ into a pair of SM fermions, only 
the decays into photons or invisible particles are possible. For mgo < Fgo the efficiency is almost 
flat. This is indicated by the dashed line which marks the average efficiency for mgo < 1 GeV. 
The shaded bands show the typical error on the efficiencies in this region. 
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Figure 4: The recoil mass spectra from =91.2 GeV a) for the decays — > e+e" and 
b) for — > OPAL data are indicated by dots with error bars (statistical error), the 

four-fermion background by the light grey histograms and the two-fermion background by the 
dark grey histograms. The dashed lines for the signal distributions are plotted on top of the 
background distributions with normalisation corresponding to the cross section excluded at 95% 
confidence level from the combination of both channels. 
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Vs=183-209 GeV, ^ e V 
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4>a (rad) lp;+P2l (GeV) 

Figure 5: Cut variables for TF* e"'"e~ at = 183-209 GeV. The OPAL data are indicated by 
dots with error bars (statistical error), the four- fermion background by the light grey histograms 
and the two-fermion background by the medium grey histograms. The signal distributions from 
a 30 GeV S° are plotted as dashed lines and those from a 100 GeV S'^ as dotted lines, respectively. 
The signal histograms are normalised to 10 and 100 times of the Standard Model Higgs-strahlung 
cross section, respectively, and show the decays S'^ gg. Each variable is shown with the cuts 
applied before the cut on this variable is done, respecting the order of cuts in Table The 
arrows indicate the accepted regions. 
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Figure 6: Cut variables for Z° ^ at ^/s 



lp^+P2l (GeV) 

183-209 GeV. The OPAL data are indicated by 



dots with error bars (statistical error), the four- fermion background by the light grey histograms 
and the two-fermion background by the medium grey histograms. The signal distributions from 
a 30 GeV S° are plotted as dashed lines and those from a 100 GeV as dotted lines, respectively. 
The signal histograms are normalised to 10 and 100 times of the Standard Model Higgs-strahlung 
cross section, respectively, and show the decays —>■ gg. Each variable is shown with the cuts 
applied before the cut on this variable is done, respecting the order of cuts in Table The 
arrows indicate the accepted regions. 
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Figure 7: The efficiency versus the mass at ^/s = 202-209 GeV for a subset of decay modes 
of a) in the Z° e"'"e~ and b) — > fi~^fj.~ channel. The minimum efficiencies which are used 
in the limits are given as well. For the other LEP 2 centre-of-mass energies the signal efficiencies 
are similar. 
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Figure 8: The recoil mass spectrum from 183-209 GeV a) for the decays Z° 



e ' e 



and b) 



for TP — > /i /i~ (lower plot). OPAL data are indicated by dots with error bars (statistical error), 
the four-fermion background by the light grey histograms and the two-fermion background by 
the medium grey histograms. The dashed lines for the signal distributions are plotted on top 
of the background distributions with normalisation corresponding to the excluded cross section 
from the combination of both channels. 
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Figure 9: Some of the cut variables for TP — > vv at ^/s = 91.2 GeV. The last plot (number of 
clusters) is one of the additional cut variables which are used in events with two electrons. The 
OPAL data are indicated by dots with error bars (statistical error), and the total background 
by the grey histograms. The distributions from a 1.2 GeV and a 30 GeV signal are plotted as 
dashed and dotted lines, respectively. The signal histograms are normalised corresponding to 
0.01 times and 0.2 times the Standard Model Higgs-strahlung cross section and show the decay 
channel — > 77. Each variable is shown with all cuts but the cut on the variable itself. The 
arrows indicate the accepted regions. 
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Figure 10: The efficiency versus the mass at ^/s = 91 GeV for the decay S°Z° — > 77Z/I7 
and S^Z" — > e'^e~i'T'. a) in Hnear mass scale and b) in logarithmic mass scale. The dashed line 
indicates the average of the efficiencies for mgo < 50 MeV. The shaded bands show the typical 
error on the efficiencies in this region. 
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Figure 11: The upper limit on the scale factor k on the cross section for the production of 
a new scalar boson in the Higgs-strahlung-process (solid line). The dot-dashed line represents 
the expected median for background-only experiments. Both limits are calculated at the 95 % 
confidence level. The dark (light) shaded bands indicate the 68% (95%) probability intervals 
centred on the median expected values. For masses mgo < 1 GeV the limits are constant. The 
lowest signal mass tested is lO^^ GeV. 
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Figure 12: Exclusion limits for the Uniform Higgs scenario at the 95% confidence level. All 
mass intervals {mA,mB) within the area bordered by the dark line are excluded from the data. 
The shaded area marks the mass points which are expected to be excluded if there were only 
background. The light grey curves indicate isolines for several values of mc- All intervals 
{rriA, m-B) to the right of each isoline are theoretically disallowed from Equation |3|. By definition, 
only intervals (mA,m-B) right to the dashed diagonal line are valid, i.e. ttt-a < ruB. 
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Figure 13: Exclusion limits on the perturbative mass scale rnc for constant K. The solid line 
represents the limits obtained from the data, and the dotted line shows the expected limit if 
there were only background. Values for mc below the lines are excluded by this analysis at the 
95 % confidence level. 
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Figure 14: Excluded parameter regions for the simplified Stealthy Higgs scenario at the 95 % con- 
fidence level. The solid line marks the region which is excluded from the data. The shaded area 
marks the region which would be excluded if the data corresponded exactly to the background- 
only prediction. The dashed lines indicate the Higgs width depending on mn and to. 
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